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Analysis  of  Self-Excited  Combustion  Instabilities  Using 

Decomposition  Techniques 


Cheng  Huang1,  William  E.  Anderson2 
Purdue  University,  West  Lafayette,  IN,  47907 

and 

Matthew  E.  Harvazinski3,  Venkateswaran  Sankaran4 
Air  Force  Research  Laboratory  (AFRL),  Edwards  AFB,  CA,  93524 


Proper  orthogonal  decomposition  (POD)  and  dynamic  mode  decomposition  (DMD)  are 
compared  with  traditional  band-pass  filtering  based  analysis  for  the  study  of  self-excited 
combustion  instabilities  in  a  longitudinal  mode  rocket  combustor.  The  POD  analysis 
approximates  the  complex  high-rank  dynamics  with  simple  lower-rank  expressions  for  the 
mode  shapes.  Each  POD  mode,  however,  is  comprised  of  multiple  acoustic  frequencies  and 
specific  modes  of  the  pressure  and  heat  release  are  not  related,  which  makes  the  analysis 
more  qualitative.  On  the  other  hand,  the  DMD  analysis  generates  a  global  frequency 
spectrum  and  each  mode  corresponds  to  a  specific  discrete  frequency.  The  DMD  result 
therefore  provide  a  quantitative  means  for  understanding  the  relationship  between  the 
pressure  modes  and  the  heat  release  modes  and  for  establishing  the  driving  mechanisms 
responsible  for  the  incidence  of  combustion  instabilities.  The  paper  uses  these  analyses  to 
describe  the  Rayleigh  index  on  a  modal  basis  to  shed  light  on  the  frequency-based  response 
of  the  combustor  flowfield. 


I.  Introduction 

Data  processing  techniques  based  on  proper  orthogonal  decomposition  (POD)  and  dynamic  model 
decomposition  (DMD)  are  powerful  and  elegant  methods  which  are  used  to  obtain  a  low-dimensional 
approximate  description  of  high-dimensional  physical  processes.  POD  has  been  extensively  employed  in  the  study 
of  turbulence  for  cold  flow  due  to  the  limitations  of  experimental  measurements  in  terms  of  both  temporal  and 
spatial  resolution.1"3  POD  can  also  be  used  to  extract  mode  shapes  and  basis  functions.  Its  application  in  combustion 
problems  has  received  some  attention  recently, 4-6  where  it  has  been  used  to  gain  insight  into  flame  dynamics.  DMD 
is  a  new  technique  that  has  not  yet  been  extensively  used  for  reacting  flow  analysis.  Thus  far  most  of  DMD’s 
applications  have  been  investigations  of  flowfield  characteristics.7'9 

Traditional  data  processing  in  combustion  instability  analysis  involves  band-pass  filtering  of  the  data  around  the 
frequencies  determined  from  a  power  spectral  density  (PSD)  analysis  of  the  raw  signal.  In  this  way,  the  correlation 
between  acoustic  and  combustion  can  be  explored  within  a  certain  frequency  range.  However,  in  many  cases, 
filtered  results  are  sensitive  to  the  filtering  methods.  Moreover,  information  about  the  couplings  between  frequencies 
might  be  lost  by  band-passing  the  signals.  One  primary  advantage  of  the  decomposition  techniques  over  filtering  is 
that  they  deal  with  the  whole  data  with  minimum  information  loss.  Moreover,  unlike  filtering  techniques,  POD  and 
DMD  do  not  require  a  prior  knowledge  or  pre-analysis  of  the  underlying  physical  and  chemical  phenomena.  Another 
motivation  is  that  they  are  capable  of  extracting  dynamically  significant  structures  from  the  flowfield  of  concern.  4-5 
Each  decomposed  mode  can  be  represented  in  terms  of  a  spatial  response  and  a  temporal  response,  which  can  in  turn 
provide  detailed  insight  into  the  dynamics  of  acoustic  and  combustion. 
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In  this  paper,  our  emphasis  is  on  the  exploratory  applications  of  POD  and  DMD  for  self-excited  combustion 
instability  in  a  gas-gas  coaxial  rocket  injector,  the  simulation  of  which  was  conducted  by  Harvzinski  et  al.10,  11 
Moreover,  distorted  triple  flame  has  been  found  using  large-eddy  simulation  of  the  same  geometry  by  Garby  et  al.12 
Experimental  measurements  of  this  configuration  have  been  published  previously  and  recent  results  are  given  by 
Feldman  et  al.13  The  present  POD  and  DMD  analyses  are  carried  out  using  the  afore-mentioned  simulation  datasets 
of  Harvzinski  et  al.10,  11  A  primary  motivating  element  for  these  studies  is  to  shed  light  on  the  fundamental 
interactions  between  the  acoustic  waves  and  the  combustion  heat  release,  which  is  a  critical  contributor  to  the 
development  of  combustion  instabilities  in  this  experiment  in  particular,  and  also  in  rocket  injector  flowfields  in 
general. 


II.  Background 

The  decomposition  analysis  generally  takes  an  ensemble  of  data,  instantaneous  or  a  series  of  snapshots,  as  inputs 
and  uses  a  matrix  processing  technique  called  singular  value  decomposition  (for  POD)  or  the  Amoldi  algorithm  (for 
DMD).  The  decomposition  of  the  original  data  sets  is  based  on  optimality  for  POD  and  frequencies  for  DMD.  Both 
techniques  result  in  modes  which  form  an  orthonormal  basis.  A  detailed  mathematical  description  of  the  POD 
technique  is  given  by  Chatterjee,14  Berkooz  et  al.,15  Cordier  and  Bergmann.16  An  introduction  to  DMD  can  be  found 
in  Rowley  et  al.17  and  Schmid.18  An  overview  of  the  mathematical  models  of  POD  and  DMD  is  given  in  the 
following  sections. 


A.  Mathematical  Model  of  Decomposition  Techniques 

The  goal  is  to  approximate  a  function,  z(x,  t ) ,  over  a  domain  of  interest  as  a  finite  sum  in  the  variables-separated 
form 

M 

k= 1 

with  the  reasonable  expectation  that  the  approximation  becomes  exact  as  M  approaches  infinity.  Note  that  in  Eq. 
(1)  there  is  no  fundamental  difference  between  x  and  t ,  but  we  usually  think  of  x  as  a  spatial  coordinate  and  t  as 
the  temporal  coordinate. 

The  representation  of  Eq.  (1)  is  not  unique.  For  example,  if  the  domain  (experimental  or  computational)  is  a 
bounded  interval  X  on  the  real  line,  then  the  functions  0^(x)can  be  chosen  as  a  Fourier  series,  Legendre 
polynomials,  Chebyshev  polynomials,  and  so  on.  For  different  selections  of  the  space-dependent  function,  O^(x), 


the  corresponding  time-dependent  function,  ak  (t) ,  will  be  different.  They  can  be  periodic  or  non-periodic,  single¬ 
frequency  dominated  or  multi-frequency  dominated. 

In  the  POD  analysis,  the  spatial  functions,  0^(x) ,  are  chosen  to  be  orthonormal  functions,  i.e. 


j  W®*,  ( x)dx 


1  if  kx  =k2 
0  otherwise 


(2) 


Then 


ak(t )  =  |  z(x,t)Ok(x)dx , 


(3) 


From  Eq.  (2)  and  (3),  it  is  noted  that,  given  that  O^(x)  is  selected  to  be  orthonormal,  the  determination  of  ak(t) 


will  be  only  dependent  on  O^(x)  rather  than  on  other  O  functions. 


B.  Singular  Value  Decomposition  (SVD)  in  POD  Analysis 

In  practice,  whole  datasets  or  snapshots  are  arranged  into  a  POD  matrix  first,  for  example,  with  each  row  containing 
the  temporal  data  and  each  column  containing  the  spatial  data.  Thus,  if  there  are  N  rows  of  temporal  data  and  m 
columns  of  spatial  data,  the  POD  matrix  will  be  of  size  Nxm  . 

Once  we  obtain  the  POD  matrix  A  ,  the  singular  value  decomposition  (SVD)  can  be  defined  as, 

A  =  UTVt  (4) 
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where  U  is  an  NxN  orthogonal  matrix,  V  is  an  m  x  m  orthogonal  matrix,  and  E  is  an  N x  m  matrix  with  all 
elements  zero  except  along  the  diagonal.  The  diagonal  elements  E..  consist  of  r  =  min(A,  m)  nonnegative  numbers 

CF. ,  defined  as  the  singular  values  of  A  .  They  are  arranged  in  decreasing  order  and  are  unique. 

In  Eq.  (4),  let  C/E  =  Q  .  Then,  A  =  QVT  .  Letting  qk  be  the  kth  column  of  Q  and  vk  be  the  kth  column  of  V  ,  we 
write  out  the  matrix  product  as 

A  =  QVT=fjqkvTk  (5) 

k= 1 

Eq.  (5)  is  the  discrete  form  of  Eq.  (1).  And  we  can  obtain  the  matrix  for  each  POD  mode, 

A  =  <iA  (6) 

In  Eq.  (6),  Ak  is  defined  as  the  kth  POD  mode  for  matrix  A  and  it  is  noted  that  Ak  has  the  same  dimensions  as  the 
matrix  A.  Moreover,  the  column  matrix  qk  can  be  interpreted  as  the  temporal  mode  of  the  kth  POD  mode,  while  the 
column  matrix  vk  represents  the  spatial  mode.  Here,  we  note  that  if  the  temporal  mode  qk  is  obtained  such  that 
C/E  =  Q  ,  the  spatial  mode  vk  would  contain  all  normalized  numbers  and  vice  versa.  Or,  we  can  write  Eq.(6)  in  this 
way, 

A  =  uk^A  (7) 

In  Eq.(7),  uk  is  the  kth  column  of  U  and  crk  corresponds  to  the  kth  singular  value  of  matrix  A  ,  which  can  further  be 
defined  as  the  mode  power  or  energy  of  the  kth  POD  mode. 

C.  Arnoldi  Algorithm  in  DMD  Analysis 

In  the  case  of  the  DMD  analysis,  in  order  to  obtain  single  frequency  dynamic  modes,  linear  mapping  is  assumed 
from  one  snapshot  to  another.  Suppose  the  data  set  is  represented  as  a  snapshot  sequence, 

K  ={V1>V2>V3> . >VA r}  (8) 

where  v.  stands  for  the  zth  flow  field  and  v.+1  =  B\i .  The  matrix  B  here  represents  the  linear  mapping  matrix. 
Therefore, 

VXN={  v19Bv19B\9 . ,Bn-\}  (9) 

Another  assumption  is  that  there  exists  a  specific  number  N ,  beyond  which  the  vector  can  be  expressed  as  linear 
combination  of  the  previous  vectors, 

=alvx+a2v2+...  +  aN_xvN_l  or  vx=VxN~la  +  r  (10) 

Hence, 

BVxn~ 1  =  V2n  =  VxN~lS  +  reN_x  (1 1) 

where, 


r0 

1  0 


5  = 


a 


\ 


10a 


(12) 


N-2 
1  aN-l 

Applying  the  eigenvalue  decomposition  for  matrix  S, 

S  =  TlAT  (13) 

where  matrix  T  is  the  eigenvector  matrix  of  S.  Suppose  a  sufficient  number  of  snapshots  are  used,  the  eigenvalues  of 
S  can  well  represent  the  eigenvalues  of  A,  which  contains  the  time-evolution  information  of  the  flowfield. 

Similarly  the  dynamic  modes  corresponding  to  single  frequency  response  can  be  constructed  as, 

*p  zz  vxN~lT-x  (14) 

Similarly,  the  original  data  set  can  be  decomposed  into  the  form  in  Eq.  (1), 


Vx  =  ^ T 


(15) 


3 

American  Institute  of  Aeronautics  and  Astronautics 


where  matrix  *¥  contains  the  dynamic  spatial  information  and  matrix  T  contains  the  temporal  evolutional 
information. 

D.  Test  Problem 

The  detailed  description  of  the  geometry  used  in  the  present  study  has  been  given  by  Harvzinski  et  al.10, 11  The 
simulation  is  concerned  with  the  combustion  instability  in  a  continuously  variable  resonant  chamber  (CVRC).  Here, 
the  geometry  used  for  the  decomposition  analysis  is  briefly  introduced  in  Fig.  1.  Decomposed  hydrogen  peroxide,  as 
the  oxidizer,  enters  from  the  left  through  the  center  oxidizer  post.  The  fuel,  gaseous  methane  enters  just  before  the 
backstep  connecting  the  post  and  the  combustion  chamber.  When  applying  the  decomposition  analysis,  information 
from  the  entire  computational  domain  is  used  for  post-processing.  The  domain  can  be  separated  into  three  regions: 
region  I  represents  the  oxidizer  post  where  there  is  only  cold  flow,  region  II  is  the  heat  release  zone,  and  region  III  is 
the  combustion  chamber.  In  the  original  simulation,  the  two-dimensional  grid  contained  200,000  grid  points  and  the 
three-dimensional  grid  contained  5,000,000  grid  points.  Information  is  interpolated  from  the  original  fine  grid  to  a 
coarser  one  to  perform  the  decomposition  analysis  due  to  the  large  memory  requirement  of  the  decomposition 
techniques  and  the  fact  that  the  resolved  boundary  layer  need  not  be  included.  Mean  values  of  pressure  and  heat 
release  are  extracted  before  post-processing.  Only  the  fluctuating  portions  of  the  pressure  and  heat  release  p  ’  and  q  \ 
are  used  and  focus  has  been  put  on  the  results  from  three-dimensional  simulation,  which  show  good  agreement  with 
the  experimental  results.13 


II 

hi 

il 

P 
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Figure  1.  Test  section  overview  including  oxidizer  post  (I),  heat  release  zone  in  the  combustion  chamber  (II) 

and  full  combustion  chamber  (III). 

Figure  2  shows  the  power  spectral  density  (PSD)  analysis  of  the  raw  pressure  signals  from  the  computational  and 
experimental  results.10  The  first  six  strongest  acoustic  frequencies  can  be  observed  from  the  figure.  According  to  the 
computational  results,  the  frequencies  corresponding  to  the  first  six  longitudinal  pressure  modes  are  1540Hz  (1L), 
3090Hz  (2L),  4636Hz  (3L),  6182Hz  (4L),  7692Hz  (5L)  and  9240Hz  (6L)  respectively.  The  experimental  results 
indicate  similar  but  slightly  lower  values.13 


3D  PSD 


Experiment  PSD 


4  6  8 

Frequency,  kHz 

Figure  2.  PSD  Analysis  for  the  raw  pressure  data  showing  the  computational  (left)  and  experimental  results 

(right)  from  Harvazinski.10 
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Figure  3  shows  the  snapshots  of  some  of  the  key  flowfield  quantities  for  one  complete  cycle,  defined  by  the  1L 
frequency  (1540Hz).  Pressure  contours  along  with  the  streamlines  are  shown  on  the  top  of  each  time  instance.  Heat 
release  contour  is  shown  on  the  bottom  left  with  spatial  mixture  fraction  and  equivalence  ratio  contours  shown  on 
the  bottom  right.  In  the  latter  case,  the  fuel  lean  region  (i.e.,  an  equivalence  ratio  less  than  unity)  has  been  blanked 
out  in  the  contour  plot.  At  the  start  of  the  cycle  (i.e.,  0°  in  Fig.  3),  a  strong  pressure  wave  is  observed  coming  from 
the  oxidizer  post  to  the  backstep  and  two  vortices  are  generated  on  the  top  and  bottom  of  the  post.  As  the  vortices 
move  downstream  ,  they  bring  along  fuel  that  mixes  with  the  oxidizer.  At  the  mid-point  of  the  cycle  (180°  in  Fig.  3), 
the  vortices  are  further  strengthened  by  the  pressure  wave  coming  upstream  in  the  chamber,  and  the  mixing  of  fuel 
and  oxidizer  is  accelerated,  which  leads  to  strong  heat  release.  As  the  pressure  wave  moves  into  the  oxidizer  post, 
the  magnitude  of  the  vortices  decreases  and  so  does  the  combustion  heat  release.  Clearly  the  acoustics  and 
combustion  are  well  correlated  with  the  significant  heat  release  occurring  when  the  pressure  is  high  at  the  combustor 
head-end. 


■°03  6  ofe  004  006 
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Figure  3.  Instantaneous  results  of  CVRC  over  one  cycle  (based  on  first  longitudinal  frequency). 


III.  Results 

The  experimental-computational  comparisons  of  the  CVRC  data  by  Feldman  et  al.13  show  that  the  three- 
dimensional  simulations  are  a  better  match  to  the  experiments  than  the  two-dimensional  axisymmetric  simulations, 
which  under-predict  the  amplitude  of  the  pressure  oscillations.  Detailed  results  using  traditional  filtering  analysis  are 
available  in  Harvazinski  et  al.  10  and  the  following  section  provides  a  brief  overview.  The  rest  of  the  discussion 
focuses  on  the  use  of  POD  and  DMD  techniques  to  shed  light  on  the  3D  simulation  results. 

A.  Results  using  Traditional  Filtering  Techniques 

Acoustic  and  heat  release  mode  shapes  are  computed  by  dividing  the  computational  domain  into  vertical  strips 
(or  discs  in  the  three-dimensional  case).  Pressure  and  heat  release  are  volume  averaged  in  each  strip.  At  the 
conclusion  of  the  simulation,  the  data  stored  in  each  strip  are  filtered  at  the  frequencies  identified  in  the  PSD  plot 
and  then  sampled  at  a  fixed  interval.  The  sampled  data  are  then  in  space  plotted  to  show  the  mode  shape  for  each 
filtered  frequency.10 
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Figure  4  shows  the  pressure  mode  shapes  for  1L,  2L  and  3L  frequencies  and  the  results  show  the  classical  mode 
shapes  in  the  chamber,  while  the  mode  shape  in  the  oxidizer  post  appears  to  be  accommodate  the  choked  condition 
at  the  inlet  as  well  as  the  area  transition  that  occurs  at  the  back-step. 


-0.10  0.00  0.10  0.20  0.30  u  w  W  w  -0.10  0.00  0.10  0.20  0.30 

X,  m  x-  m  X,  m 

Figure  4.  Pressure  mode  shapes  obtained  by  filtering  around  1L  (left),  2L  (middle)  and  3L  (right) 

frequencies  from  Harvazinski. 10 


Figure  5  shows  the  heat  release  mode  shapes.  The  first  three  modes  are  dominant  and  show  that  the  peak  heat 
release  is  concentrated  near  0.04m  downstream  from  the  backstep.  For  the  1L  mode  a  single  peak  is  observed.  For 
the2L  mode  an  additional  small  peak  appears  just  downstream  of  the  backstep,  while  for  the  3L,  yet  another  peak 
appears  downstream  of  the  dominant  peak. 


0.10 
x,  m 


0.10 

x,  m 


0.05  0.10 

x,  m 


Figure  5.  Heat  release  mode  shapes  obtained  by  filtering  around  1L  (left),  2L  (middle)  and  3L  (right) 

frequencies  from  Harvazinski.10 


B.  POD  Results 

Figure  6  shows  the  POD  mode  power  as  a  function  of  the  mode  number.  We  note  that  the  mode  power  is  given 
by  the  value  of  ok  in  Eq.(7),  which  indicates  the  significance  of  each  POD  mode.  The  POD  mode  power  distribution 
of  p'  has  a  sharp  decrease  after  the  first  mode  and  decays  to  zero  very  quickly  as  the  mode  number  increases. 
However,  the  power  distribution  of  the  q'  mode  is  different  in  that  it  decreases  much  more  gradually.  The  first  POD 
mode  of  p'  contains  11%  of  the  total  information  while  that  of  q'  has  only  0.66%.  In  the  figure,  the  summation  of 
first  200  POD  modes  of  p'  and  q  ’  are  given  by  the  solid  line.  For  p\  the  first  200  modes  account  for  approximately 
55%  of  the  energy,  while  for  q\  the  first  200  modes  account  for  only  about  22%  of  the  energy.  This  indicates  that 
more  POD  modes  of  q'  are  need  to  fully  represent  the  behavior  of  the  heat  release  rate.  The  source  of  this  behavior 
can  be  examined  by  looking  at  the  PSD  plot,  shown  in  Fig.  7.  In  particular,  the  PSD  of  the  heat  release  perturbations 
shows  lower  power  at  the  dominant  frequencies,  which  may  be  attributed  to  the  richer  characteristics  of  the 
combustion  process  compared  with  the  acoustic  response  which  is  more  predictable  and  organized. 
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Figure  6.  POD  mode  power  and  accumulated  energy  of  pressure  oscillations  (left)  and  heat  release 

oscillations  (right). 
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Figure  7.  Power  spectral  density  of  pressure  (top)  and  heat  release  (bottom)  oscillations  near  the  location 

of  maximum  heat  release. 


In  the  next  two  sections,  mode  shapes  of  the  pressure  and  heat  release  oscillations  as  given  by  the  POD  analysis 
are  discussed.  Both  the  spatial  contours  and  temporal  modes  are  shown  for  each  POD  mode.  Since  only  the 
fluctuation  terms  are  included  in  the  decomposition  analysis,  negative  values  shown  in  the  figures  below  indicate  a 
decrease  in  magnitude  relative  to  the  mean  and  vice  versa. 
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a.  POD  Modes  of  Pressure  Oscillations 

The  first  six  POD  modes  of  the  pressure  oscillations  for  the  combined  combustion  chamber  and  oxidizer  post 
configuration  are  shown  in  Figs.  8-10.  In  each  case,  we  show  the  spatial  modes  (top)  as  well  as  the  temporal  mode 
and  the  PSD  of  the  temporal  mode  (bottom).  Contours  of  the  POD  spatial  modes  are  plotted  on  a  slice  which  goes 
through  the  center  of  the  domain.  Figure  8  shows  that  the  dominant  frequency  is  1550Hz  which  corresponds  to  the 
first  longitudinal  mode.  The  first  POD  mode  is  also  subject  to  minor  influences  from  the  2L,  3L,  4L,  and  5L 
acoustic  modes  as  evident  in  the  PSD  plot  of  temporal  mode. 

Multiple  frequencies  continue  to  appear  in  the  higher  POD  modes.  The  higher  POD  modes  tend  to  be  dominated 
by  the  higher  acoustic  frequencies  present  in  the  chamber.  This  is  consistent  with  the  energy  distribution  plot  shown 
in  Fig.  6.  Since  the  POD  is  designed  to  obtain  an  optimal  representation  of  the  original  data  set,  it  can  be  concluded 
that  the  proper  description  of  the  acoustic  responses  within  the  computational  geometry  is  complex  and  cannot  be 
approximated  using  a  single  frequency.  Instead,  a  combination  of  multiple  frequencies  is  required.  Even  though  the 
same  frequencies  appear  in  multiple  temporal  POD  modes,  their  relative  contributions  are  different  for  each  mode 
because  of  the  different  spatial  modes  present.  For  instance,  the  3rd  and  4th  POD  modes  in  Fig.  9  is  dominated  by  the 
2L  and  3L  modes,  while  the  5th  POD  mode  in  Fig.  10  is  dominated  by  the  3L,  4L  and  5L  modes  and  the  6th  POD 
mode  (also  in  Fig.  10)  is  dominated  by  the  5L  and  6L  modes. 
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Figure  8.  POD  analysis  of  pressure  oscillations  from  3D  simulations  with  1st  mode  (left),  2nd  mode  (right). 
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Figure  9.  POD  analysis  of  pressure  oscillations  from  3D  simulations  with  3rd  mode  (left),  4th  mode  (right). 
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Figure  10.  POD  analysis  of  pressure  oscillations  from  3D  simulations  with  5th  mode  (left),  6th  mode  (right). 


b.  POD  Modes  of  Heat  Release  Oscillations 

The  first  six  POD  modes  of  heat  release  oscillations  are  shown  in  Figs.  11  -  13.  Again,  in  each  case,  we  show  the 
spatial  mode  as  well  as  temporal  mode  and  the  PSD  of  the  temporal  mode.  Most  of  the  combustion  occurs  in  section 
II  (i.e.,  the  heat  release  zone  in  Fig.  1)  and,  for  this  reason,  the  heat  release  POD  modes  are  shown  only  for  this 
section.  The  first  POD  mode  of  the  heat-release  in  Fig.  11  shows  a  combustion  pulse  between  0.03  and  0.04m.  The 
second  heat-release  POD  mode  shows  a  half-wave  configuration  in  terms  of  the  spatial  mode  contour.  As  in  the  case 
of  pressure,  each  POD  mode  is,  in  fact,  a  combination  of  longitudinal  mode  frequencies.  The  frequency  values 
correspond  to  the  pressure  modes  very  well;  however,  the  specific  combination  of  longitudinal  modes  in  each  case  is 
different  from  the  pressure  case. 

It  is  interesting  that  the  combustion  heat  release  in  the  1st  POD  mode  in  Fig.  11  appears  to  come  from  the  main 
vortex  in  Fig.  3,  which  we  noted  earlier  is  strengthened  by  the  incoming  pressure  wave.  This  vortex  accelerates  the 
mixing  of  fuel  and  oxidizer  and  promotes  the  heat  release.  Moreover,  the  location  where  the  vortex  meets  the  strong 
pressure  wave  is  consistent  with  the  location  of  the  heat-release  spatial  mode.  This  result  is  therefore  a  clear 
indication  of  the  combustion-pressure  coupling  that  plays  a  critical  role  in  generating  and  sustaining  the  combustion 
instabilities.  As  the  vortex  moves  downstream,  it  starts  to  be  influenced  by  the  additional  acoustic  modes,  which 
separates  the  single  main  vortex  into  several  smaller  ones.  This  can  be  observed  from  the  2nd  POD  mode  also  shown 
in  Fig.  11.  Here,  an  additional  combustion  heat-release  zone  appears  further  downstream,  which  is  opposite  in  phase 
with  the  primary  heat-release  location.  As  shown  in  the  PSD  plot  in  1st  and  2nd  POD  modes  in  Fig.  11,  these  two 
combustion  modes  are  highly  responsive  to  the  first  three  acoustic  modes.  Thus,  the  two  combustion  modes  appear 
to  represent  the  major  combustion  dynamics  fluctuations  that  we  observed  earlier  in  Fig.  3. 
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Figure  11.  POD  analysis  of  heat  release  oscillations  from  3D  simulations  with  1st  mode  (left),  2nd  mode 

(right). 


After  the  high  amplitude  pressure  wave  passes  through  the  chamber  and  enter  the  post,  higher  acoustic  modes 
with  lower  magnitudes  come  into  play,  perhaps  influenced  by  the  dynamics  in  the  post.  These  higher  modes  lead  to 
additional  vortex  splitting  and  in  turn  lead  to  lower  amplitude  combustion  fluctuations,  which  are  captured  by  the  3rd 
and  4th  POD  modes  in  Fig.  12.  These  two  modes  describe  the  downstream  transition  of  the  smaller-scale  combustion 
and,  moreover,  they  suggest  an  extra  upstream  combustion  location.  In  Fig.  3  (at  240°  and  300°),  small  vortices  are 
observed  to  be  pushed  upstream  towards  the  backstep  by  the  upcoming  pressure  waves.  This  results  in  enhanced 
mixing  of  fuel  and  oxidizer  near  the  backstep  location.  Moreover,  the  PSD  plots  indicate  that  the  4L,  5L  and  6L 
modes  are  dominant  in  these  combustion  responses.  Further  description  of  the  role  of  longitudinal  frequencies  in  the 
combustion  dynamics  is  discussed  later  in  section  C  as  part  of  the  DMD  analysis. 
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Figure  12.  POD  analysis  of  heat  release  oscillations  from  3D  simulations  with  3rd  mode  (left),  4th  mode 

(right). 


The  5th  and  6th  POD  mode  of  the  heat  release  in  Fig.  13  do  not  suggest  a  strong  repeatable  behavior  in  terms  of 
the  temporal  modes  shown  in  the  bottom  left  although  they  show  some  peaks  around  1L,  2L  and  3L  frequencies  in 
terms  of  the  PSD  plots  in  bottom  right.  Moreover,  the  spatial  modes  do  not  indicate  a  predictable  pattern  as  do  the 
first  four  POD  modes.  In  fact,  there  is  some  evidence  of  radial  or  azimuthal  modes  in  these  patterns  which  should  be 
investigated  further. 
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Figure  13.  POD  analysis  of  heat  release  oscillations  from  3D  simulations  with  3rd  mode  (left),  4m  mode  (right). 
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Although  in  terms  of  the  POD  mode  power  spectrum  in  Fig.  6,  the  heat  release  oscillations  require  a  large 
number  of  POD  modes,  we  note  that  the  primary  combustion  response  is  captured  by  the  first  few  POD  modes.  The 
higher  POD  heat  modes  can  be  regarded  as  representations  of  spontaneous  and  chaotic  combustion  processes  which 
are  not  necessarily  coupled  with  the  acoustic  modes.  This  indicates  that  focusing  on  the  lower  POD  heat  release 
modes  is  sufficient  to  shed  insight  into  the  physical  phenomena  of  relevance  to  combustion  stability. 


C.  DMD  Results 

Unlike  the  POD  analysis,  the  DMD  technique  provides  modes  at  discrete  frequencies.  Figure  14  shows  the  DMD 
mode  power  versus  frequencies  allowing  the  importance  of  each  frequency  to  be  determined.  Both  pressure  and 
heat-release  modes  are  shown.  As  expected,  the  heat  release  oscillations  show  strong  response  to  the  longitudinal 
acoustic  frequencies.  But,  other  than  the  expected  acoustic  frequencies,  there  are  also  additional  peaks  of  lower 
magnitude  in  the  high  frequency  region  from  10,000  to  20,000  Hz.  This  confirms  what  has  been  concluded 
previously  from  the  POD  mode  results,  that  in  addition  to  the  pressure  mode  response,  the  combustion  dynamics 
also  contains  additional  physics  that  can  be  characterized  as  the  result  of  the  turbulent  reacting  flowfield. 
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Figure  14.  DMD  mode  power. 
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Mode  studies  of  the  pressure  and  heat  release  oscillations  using  DMD  analysis  are  shown  in  Figs.  15  -  17.We 
note  that  the  DMD  temporal  modes  are  very  simple  since  they  include  only  a  single  frequency.  For  this  reason,  only 
spatial  modes  are  shown.  The  pressure  modes  produced  by  the  DMD  are  immediately  recognizable  and  look  like  the 
acoustic  modes  one  would  expect  in  the  chamber.  Figure  15  shows  the  DMD  pressure  modes  for  the  first  two 
longitudinal  modes.  In  the  combustor,  the  first  mode  is  a  half-wave  and  the  second  mode  is  a  full-wave.  The  mode 
shapes  within  the  oxidizer  plot  are  slightly  more  complicated  but  also  correspond  to  the  above  trend.  Specifically, 
the  combined  mode  shapes  also  correspond  to  the  band-passed  one-dimensional  mode  shapes  derived  in  Fig.  4. 

DMD  modes  of  the  heat  release  oscillations  in  response  to  1L  and  2L  acoustic  frequencies  in  Fig.  15  appear  to 
give  similar  interpretations  as  in  the  case  of  the  POD  results  given  in  Fig.  11.  Importantly,  the  1L  DMD  heat-release 
mode  shows  a  pulsating  response  near  0.04m,  which  confirms  that  this  response  is  mainly  driven  by  1L  acoustic 
frequency.  The  2L  DMD  heat  release  mode  demonstrates  that  the  combustion  starts  to  function  in  alternate 
locations,  both  upstream  and  downstream  of  the  primary  response  region  of  the  first  mode.  We  note  that,  in  contrast, 
the  POD  modes  represent  the  response  to  a  combination  of  acoustic  modes  and,  therefore,  it  is  difficult  to  discern  the 
precise  driving  mechanisms.  Thus,  the  DMD  analysis  provides  more  precise  mode-by-mode  representation  of  the 
driving  mechanisms  and  is  therefore  a  more  quantitative  approach  than  the  POD  analysis. 


Figure  15.  DMD  modes  of  pressure  and  heat  release  with  1L  frequency  (left)  and  2L  frequency  (right). 


When  it  comes  to  responses  to  higher  acoustic  frequencies  shown  in  Fig.  16  and  Fig.  17,  the  DMD  analysis  gives 
the  expected  pressure  mode  shapes  in  the  combustion  chamber  (one  and  half  waves  for  3L,  two  full  waves  for  4L, 
two  and  half  waves  for  5L  and  three  full  waves  for  6L)  while  in  the  oxidizer  post,  there  appear  to  be  some  evidence 
of  mode-shifting  with  the  modes  in  the  chamber.  For  example,  the  DMD  pressure  mode  of  the  3L  frequency  shows  a 
full  wave  in  the  oxidizer  post.  This  lack  of  tuning  might  be  attributed  to  the  different  sound  speeds  in  the  chamber 
and  the  oxidizer  post  because  of  the  lower  temperature  of  the  unbumed  propellants  in  the  post. 

The  3L  DMD  heat-release  mode  in  Fig.  16  (left)  suggests  a  similar  spatial  configuration  as  the  2L  DMD  heat- 
release  mode.  This  indicates  that  the  2L  and  3L  frequencies  contribute  significantly  to  the  alterations  of  the 
combustion  processes  between  0  m  and  0.08m  in  the  chamber  while  the  response  driven  by  the  1L  frequency  is 
centered  at  0.04m.  The  higher  DMD  heat  release  modes  in  Fig.  16  (right),  Fig.  17  (left)  and  Fig.  8  (right)  show 
similar  results  although  the  heat  release  zones  are  observed  to  be  more  compact  presumably  because  of  the  higher 
frequencies  (and,  consequently  shorter  wavelengths)  involved. 
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Figure  16.  DMD  modes  of  pressure  and  heat  release  with  3L  frequency  (left)  and  4L  frequency  (right). 


Figure  17.  DMD  modes  of  pressure  and  heat  release  with  5L  frequency  (left)  and  6L  frequency  (right). 


Other  than  the  longitudinal  acoustic  frequencies  discussed  so  far,  the  DMD  analysis  also  indicates  three  lower 
amplitude  peaks  between  the  longitudinal  frequencies,  highlighted  in  the  spectrum  in  Figure  14  at  14004,  14371  and 
14538  Hz  respectively.  Spatial  contours  in  response  to  these  three  frequencies  are  shown  in  Figure  18.  In  terms  of 
the  pressure  contours,  14004Hz  corresponds  to  the  1st  tangential  (IT)  acoustic  mode,  14371Hz  appears  to  show  the 
1T2L  mode,  and  14538Hz  appears  to  correspond  to  the  1T3L  mode.  These  tangential  modes  can  be  theoretically 
described  by  the  linear  wave  equation  but  have  not  been  previously  observed  in  the  CVRC  simulations  using 
traditional  data-processing  techniques.  Heat  release  in  response  to  these  high  frequencies  is  dominated  by  small 
scale  combustion,  which  may  be  more  related  to  turbulence  than  acoustics.  In  fact,  the  DMD  spectrum  displays  no 
strong  peaks  for  the  heat  release  at  these  frequencies,  which  allows  us  to  conclude  that  the  tangential  acoustic  modes 
play  weaker  role  in  influencing  combustion  than  do  the  longitudinal  modes.  Nevertheless,  these  results  demonstrate 
the  capability  of  the  DMD  analysis  to  shed  light  on  underlying  physical  phenomena  that  are  not  otherwise  evident. 


13 

American  Institute  of  Aeronautics  and  Astronautics 


Figure  18.  DMD  modes  of  pressure  and  heat  release  with  14004Hz  (top  left),  14371Hz  (top  right)  and 

14538  Hz  (bottom). 

At  this  point,  we  can  also  draw  some  conclusions  on  the  relative  utility  of  the  POD  and  DMD  analysis.  The  POD 
technique  is  designed  to  extract  the  representative  dynamic  modes  out  of  the  original  dataset  or  snapshots.  It  offers  a 
mathematically  optimal  lower-order  approximation  of  the  original  high-order  system  and  therefore  may  provide  a 
good  quantitative  understanding  of  the  observed  phenomena.  Connections  between  frequencies  can  be  studied 
because  in  most  cases,  the  POD  modes  are  multi-frequency  driven.  However,  the  couplings  between  different 
variables,  such  as  pressure  and  heat  release,  are  not  easily  determined.  Specifically,  the  POD  modes  for  pressure  and 
heat  release  do  not  have  a  one-to-one  correspondence  and  may  in  general  contain  the  responses  to  different 
underlying  mode  structures.  On  the  other  hand,  the  DMD  technique  provides  more  quantitative  interpretations  of  the 
response  to  a  specific  frequency.  Pressure  and  heat  release  can  be  coupled  easily  on  the  basis  of  frequencies,  which 
can  be  helpful  in  developing  fundamental  insight  into  the  coupling  mechanisms  responsible  for  driving  the 
instabilities.  Moreover,  the  approach  can  be  useful  for  developing  appropriate  transfer  functions  between  the 
acoustics  and  the  combustion,  which  can  be  used  in  lower-order  engineering  models  of  combustion  instability. 

D.  Applications  of  Decomposition  Techniques  in  Calculating  Rayleigh  Index 

The  Rayleigh  index  is  defined  as  the  degree  of  correlation  between  the  pressure  and  the  combustion  heat  release. 
The  Rayleigh  criterion  states  that  when  the  pressure  and  heat  release  are  in-phase  (i.e.,  a  positive  Rayleigh  index), 
combustion  instabilities  are  driven,  while  when  they  are  out-of-phase  (i.e.,  a  negative  Rayleigh  index),  combustion 
instabilities  are  damped.  The  Rayleigh  index  is  therefore  a  useful  diagnostic  tool  for  obtaining  a  quantitative 
understanding  of  combustion  instability  incidence.  Previous  studies  of  the  Rayleigh  index  are  based  on  filtered 
signals  around  the  concerned  acoustic  frequencies  as  reported  by  Harvazinski  et  al.19  Sample  three-dimensional 
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results  of  the  Rayleigh  index  using  filtered  data  for  the  first  three  acoustic  frequencies  have  been  analyzed  by 
Harvazinski10  and  are  reproduced  in  Fig.  19.  We  will  discuss  these  trends  in  the  context  of  similar  results  obtained 
using  the  POD  and  DMD  analysis. 


Figure  19.  Rayleigh  index  calculated  from  filtering  techniques  by  Harvazinski. 10 


The  Rayleigh  index  calculated  from  the  POD  modes  is  shown  in  Fig.  20.  Because  the  decomposition  analysis  is 
performed  independently  for  each  variable  (p  ’  and  q  %  it  is  hard  to  correlate  pressure  and  heat  release  oscillations 
using  the  POD  results.  As  discussed  above,  each  POD  mode  does  not  necessarily  correspond  to  a  single  frequency. 
Therefore,  to  calculate  the  Rayleigh  index  from  the  POD  modes,  the  non-filtered  raw  signal  of  heat  release 
oscillations  is  used  rather  than  the  filtered  or  decomposed  ones.  Moreover,  the  decomposed  pressure  POD  modes  are 
used  in  terms  of  the  dominant  frequencies.  Based  on  the  results  in  Figs.  8-10,  the  1st  ,3rd  and  5th  POD  pressure 
modes  have  dominant  responses  from  the  1L,  2L  and  3L  acoustic  frequencies  respectively.  For  this  reason,  we  select 
these  three  POD  modes  to  derive  the  Rayleigh  index.  In  Fig.  20,  the  contour  scales  are  normalized  by  the  strongest 
Rayleigh  index  magnitude,  which  comes  from  1st  POD  mode  of  pressure  in  these  three  cases. 

Comparing  the  results  in  Fig.  19  and  Fig.  20,  the  Rayleigh  index  due  to  the  1st  POD  pressure  mode  (which  is 
mainly  driven  by  1L  frequency)  gives  similar  interpretations  of  the  combustion  heat-release  mechanism.  A  dominant 
driving  region  lies  downstream  near  the  0.04  m  location  with  a  small  damping  region  sitting  close  to  0.02m.  The 
corresponding  results  for  the  3  rd  POD  mode  shows  that  the  maximum  Rayleigh  index  is  only  10%  compared  to  the 
first  mode  result;  however,  it  shows  the  same  locations  for  the  strongest  driving  and  damping  regions  as  do  the 
filtered  results  in  Fig.  19.  Finally,  the  Rayleigh  index  results  for  the  5th  POD  pressure  mode  appear  to  be  similar  to 
the  3L  filtered  Rayleigh  index  with  a  relatively  low  damping  region  near  the  inlet  (backstep)  and  a  strong  driving 
area  0.02m  downstream. 

Unlike  the  POD  modes,  the  DMD  modes  of  the  pressure  and  heat  release  oscillations  can  be  correlated  by 
frequency.  Figure  21  shows  the  Rayleigh  indices  based  on  the  DMD  modes  corresponding  to  the  1L,  2L  and  3L 
acoustic  frequencies.  Comparing  the  results  shown  in  Fig.  19,  Fig.  20  and  Fig.  21,  the  Rayleigh  index  driven  by  1L 
frequency  predicts  a  single  dominant  driving  location  centered  around  0.04  m,  except  that  in  the  filtered  results,  the 
location  is  slightly  downstream  in  the  filtered  case.  The  2L  and  3L  DMD  Rayleigh  indices  are  also  qualitatively 
similar  to  the  results  observed  in  the  filtering  and  POD  analysis  in  that  they  predict  approximately  the  same 
locations  for  driving  and  damping  mechanisms. 
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Figure  20.  Rayleigh  index  calculated  from  POD  modes. 


Figure  21.  Rayleigh  index  calculated  from  DMD  modes, 
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Overall,  the  Rayleigh  index  results  given  by  the  DMD  analysis  are  qualitatively  similar  to  the  filtered  results,  but 
they  are  not  identical.  This  is  because  when  filtering  the  signals,  the  results  can  be  influenced  by  the  way  the  signals 
are  processed.  For  example,  different  bandwidths,  different  types  of  filtering  windows  or  different  orders  of  filtering 
can  influence  the  results.  Therefore,  filtering  cannot  precisely  ensure  a  purely  single-frequency  response.  In  that 
respect,  the  DMD  results  are  perhaps  the  most  reliable  way  to  quantify  such  effects. 

IV.  Conclusions 

This  paper  is  concerned  with  the  analysis  of  three-dimensional  simulation  results  for  self-excited  combustion 
instabilities  in  a  longitudinal  mode  rocket  chamber.  Both  proper  orthogonal  decomposition  (POD)  and  dynamic 
mode  decomposition  (DMD)  techniques  are  used  as  data-processing  methods.  The  overall  goal  is  to  use  these 
techniques  to  shed  light  on  the  fundamental  coupling  mechanisms  between  the  chamber  acoustics  and  the 
combustion  dynamics  and  to  thereby  understand  the  driving  mechanisms  that  promote  the  incidence  of  combustion 
instabilities. 

POD  analysis  provides  a  series  of  POD  modes  of  decreasing  energy  content.  However,  each  of  the  modes 
correspond  to  a  mix  of  acoustic  modes  and,  moreover,  there  is  no  direct  connection  between  the  individual  pressure 
and  heat  release  modes.  These  factors  make  the  POD  analysis  more  useful  as  a  means  of  lower-order  representation 
of  the  solution,  but  the  analysis  is  more  limited  in  terms  of  establishing  a  quantitative  basis  for  the  combustion- 
acoustic  coupling  mechanisms.  The  frequency-based  DMD  technique,  on  the  other  hand,  provides  modes  that 
correspond  to  single  frequencies.  Moreover,  the  corresponding  pressure  and  heat-release  modes  contain  useful 
information  regarding  how  these  two  variables  are  coupled  to  each  other  and  therefore  provide  a  more  quantitative 
basis  for  understanding  combustion  instability. 

In  terms  of  the  heat  release  fluctuations,  both  POD  and  DMD  analysis  clearly  reveal  a  dominant  pulsating 
combustion  mode,  which  is  driven  mainly  by  the  1L  acoustic  mode.  The  corresponding  Rayleigh  index  derived  from 
the  DMD  analysis  shows  that  this  combustion  response  is  in  phase  with  the  pressure  and  leads  to  a  strong  positive 
Rayleigh  index  that  is  centered  around  the  0.04  m  location  downstream  of  the  back- step.  The  2L  and  3L  DMD  heat 
release  modes  show  that  there  are  additional  combustion  response  locations  on  either  side  of  the  dominant  mode. 
The  corresponding  Rayleigh  index  reveals  both  damping  and  growth  regions,  although  the  growth  regions  are 
considerable  weaker  than  that  obtained  for  the  1L  mode.  This  is  to  be  expected  since  the  present  case  corresponds  to 
a  1L  mode  instability  in  the  chamber. 

Future  work  will  include  extending  these  studies  to  more  stable  regimes  of  combustor  operation  and  using  the 
POD  and  DMD  analyses  to  predict  the  damping  mechanisms  that  are  dominant  under  stable  operating  conditions. 
Further,  the  use  of  DMD  analysis  to  derive  combustion  response  functions  will  also  be  investigated. 
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